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A numerical analysis is made of the breakdown of the quan- 
tum Hall effect caused by the Hall electric field in competition 
with disorder. It turns out that in the regime of dense impu- 
rities, in particular, the number of localized states decreases 
exponentially with the Hall field, with its dependence on the 
magnetic and electric field summarized in a simple scaling law. 
The physical picture underlying the scaling law is clarified. 
This intra subband process, the competition of the Hall field 
with disorder, leads to critical breakdown fields of magnitude 
of a few hundred V/cm, consistent with observations, and ac- 
counts for their magnetic-field dependence oc B 3 ' 2 observed 
experimentally. Some testable consequences of the scaling law 
are discussed. 

73.40.Hm,73.20.Jc 



I. INTRODUCTION 

Since the discovery of the quantum Hall effectEftl 
(QHFJ .there has been considerable interest in its break- 
downEniij due to high currents. Rich accumulation of ex- 
perimental data has recently shed new light on this break- 
down phenomenon: A series of experiments by Kawaji et 
alM using butterfly-shaped Hall bars clearly showed that 
the breakdown of the QHE is controlled by the current 
density or the Hall field rather than the current itself 
and that for samples with a variety of electron concen- 
tration and mobility the critical values of the Hall field 
at different diagonal-resistance plateaus scale like oc B 3 / 2 
with the magnetic field B. This B 3 ! 2 dependence disfa- 
vors some of the breakdown mechanisms proposed earlier, 
such as thpse attributed to .electron heatingQ and phonon 
emission^ Zener tunnelingErtLatll of electrons across the 
Landau gap appears consistent with it but theoretically 
yields critical fields one order of magnitude larger than 
observed experimentally. 

In a previous paperLLa we presented a numerical study 
of current distributions in Hall bars. There we noted, as 
a by-product of the study, the possibility that the Hall 
field in competition with disorder in the sample interior 
would be responsible for the breakdown of the QHE, with 
a simple estimate of the competition that leads to criti- 
cal fields of magnitude consistent with observations and 
their B 3 / 2 scaling law. The purpose of the present pa- 
per is to verify by a detailed numerical analysis that this 
intrasubband process can indeed account for the break- 
down phenomenon. Actually Stormer et alia earlier noted 
that the action of a Hall field could qualitatively explain 
many of their experimental findings. Recently Shimada 



et alx3 also suggested, through a measurement of activa- 
tion energy, that the J3 3 / 2 scaling law is related to field- 
dependent broadening of the spectra of extended states. 
We shall study numerically how the fraction of delocal- 
ized states per Landau subband depends on the electric 
and magnetic field and clarify the physical picture un- 
derlying the scaling law. .— . 

In Sec. II we review a theoretical frameworkcS for our 
numerical simulations. We examine the case of a single 
impurity in Sec. Ill and handle Hall samples with random 
impurities in Sec. IV. Section V is devoted to a summary 
and discussion. 



II. FORMALISM 

Consider electrons confined to a strip of length L x and 
width L y (or formally, a strip bent into a loop of circum- 
ference L x ), described by the Hamiltonian: 

H = H + U(x,y)-eA (y), (2.1) 

Ho = lu,{fpl + (l/£ 2 )(y-yo) 2 }, (2.2) 

written in terms of u = eB/m, the magnetic length £ = 
1/veB and yo = p x /(eB). Here a uniform magnetic field 
B normal to the plane has been introduced by use of the 
Landau-gauge vector potential (— By,0). 

We shall detect the current j x (x, y) flowing in the pres- 
ence of a Hall potential Aq (y) and an impurity potential 



U{x, y) = ^2 ^ S ( x ~ x i) S (v 



Vi 



(2.3) 



which consists of short-range impurities of strength Ai 
distributed over the sample at position (Xi,yi). For sim- 
plicity, we here take Ao(y) — —yE y which produces a 
uniform Hall field E y in the y direction. 

The cyclotron motion of an electron is described by 
H , whose eigenstates are Landau levels \N) = \n,yo) 



labeled by integers n = 0,1,2,- 
wave functions of the form 



(x,y\N) = (L x ) 



-1/2 izp a 



and j/o 



4>n(y,yo)- 



'Px 



with 



(2.4) 



The 4>n(y,yo) are highly localized around y ~ yo with 
spread Ay ~ 0(£), and are given by the usual harmonic- 
oscillator wave functions for electrons in the sample bulk. 
In what follows we shall take no explicit account of the 
electron edge states, which, being always extended, have 
little to do with field-induced derealization. 



With (x,y\N) taken to be periodic in x, yo — £ 2 p x is 
labeled by integers k: 



y = (2ir£ 2 /L x )k = y [k}. 



(2.5) 



We shall henceforth use this discrete label k rather than 
yo to refer directly to each electron state, with normal- 
ization (N\N') = d nn 'S k k'. 

For numerical analyses it is advantageous to handle the 
Hamiltonian Hmn 1 = (iV|H|iV) in N — (n,k) space: 

H NN > = uj(n+ ^)S NN >+ U NN > + eE y y NN >, (2.6) 

U NN , = ±J2 X i e ~ i(y ™'° )Xi/i2 tn(yiM)ct>n>(yiW' ), (2-7) 



Vnn> 

- 1 - n. n 



{yoK 



+ lY nn ,}5, 



kk' 



y/n'/2 6 n +i, n > + y/n/2 S. 



n,n' + 1 j 



(2.8) 
(2.9) 



where y = y [k] and y' = y [k']. 

When the disorder Unn 1 is weak compared with the 
level gap, i.e., \Xi/2-K£ 2 \ < u, one can diagonalize Hmn' 
with respect to Landau-level labels (n, n') by a suitable 
unitary transformation H w — \VHW~ 1 . In particular, 
a simple 0(U) choice W 



,,iA 



with 



ih-NN' 



lEy Y nn i 5kk' + Unn 1 
u(n — n') 



(for n ^ ri), (2.10) 



leads to the Hamiltonian (N\H W \N') = S nn ,H^ n (k,k') 
governing each impurity-broadened Landau subband: 



H, 



ir 



(k, k') = [oj(n+ |) + eE y yo^Skk' 

+U nn {k,k') + 0(U 2 /Lu), 



(2.11) 



where we have set Unw —* U nn i (k, k'). 

Numerically diagonalizing i/^(fc, k'), one obtains 
eigenstates \a) forming the nth subband, with wave func- 
tions of the form (n, k\a). Then the current density car- 
ried by an eigenstate \a) is calculated through 

ji a) (x,y) = eujRe [(a\W\x){x\{y - y )W- 1 \a)] . (2.12) 

Impurities capture orbiting electrons and make them 
localized in space. In contrast, an applied (or internally 
generated) Hall field causes a drift of electrons with veloc- 
ity v x ~ Ey/B and works to delocalize them. A simple 
estimate of energy costll3 reveals how these two effects 
compete: Consider an impurity of strength A placed at 
x = y = on an infinite plane with E y = 0. It traps 
an electron and in each Landau level n there arises one 
localized state of energy shift Ae w X/(2tt£ 2 ) relative to 
the level center uj(n + ^). [This is readily seen if one notes 
that U nn (k, k') in the present case is a projection to the 
nth level (oc |n)(n|) of a harmonic oscillator with coordi- 
nate yo-] It is thus convenient to define a dimensionless 
strength s by A/(2ir£ 2 ) — suj. As a Hall field is turned 
on, this localized state will perceive over its spatial ex- 
tent of 0(£) an energy variation of magnitude ~ e£E y . 
Accordingly, if the field becomes so strong that 



ei^l >|A|/(27r^ 2 ) = | s | w , 



(2.13) 



the impurity would fail to capture any electron. In the 
next section we shall make this criterion for derealization 
quantitative by a numerical analysis . 



III. THE CASE OF A SINGLE IMPURITY 



Consider a sample of length L x = \ / 2tt£ x 16 w 40£ 
and width L y = \/2tt£ x8« 2(W, supporting 16 x 8 = 
128 electron states per level with —\L y < y < \L yi 
or y = (2ir£ 2 /L x )k with k = -64, -63, • • • ,63. Place 
an impurity of strength s = 0.1 at (x, y) — (W£, 0) in 
the middle of the sample and examine the electron states 
about the impurity by varying the field E y or the ratio 



R = e£\E y \/suj. 



(3-1) 



It turns out that the electrons states near the boundary 
y = ±\L y ~ ±10^ are hardly affected by the impurity. 

A close look into the n = 0,1,2,3 levels reveals the 
following general features: In the nth level there arise 
(n + 1) states of sizable spread (in yo) ~ 0(£), all of 
which turn out localized when E y is sufficiently weak 
(e.g., R < 0.001). Only one of them, having Ae w suj, is 
tightly localized. The other n states gain practically no 
energy shift Ae rs and get easily delocalized as E y is 
increased, in accordance with the criterion (2.13). 

In what follows we shall focus on the tightly-localized 
states with Ae w suj; there is one such state per level. 
They carry no current while E y is sufficiently weak. For 
R > 0.15 they start carrying an appreciable amount of 
current ,P° C = J dy j l ° c (x,y), increasing with E y on av- 
erage and becoming sizable at some discrete values of 
E y . Figure 1 shows such current-field characteristics for 
a potentially-localized state in the n = level, along with 
a smooth curve interpolating through the local minima. 

A direct look into current distributions shows that for 
these discrete values of E y a localized state strongly inter- 
feres with an extended state nearly degenerate in energy 
and lying within its spatial extent. Noting that extended 
states are equally spaced with separation Ayo = 2n£ 2 /L x 
in t/o space, one expects that such interference takes place 
when suj w eEy/Sy^y. integer, i.e., at discrete values of 
the ratio R, 



1?" 



(£/Ay )/K 



(3.2) 



with some integers K = 1, 2, • • •. Indeed, the sequence 

R disc w 0.1705, 0.1803, 0.1913, 0.204, • • • , (3.3) 

obtained for the localized state of the n = level is 
reproduced from Eq. ( fL2; ) with Ay rs 0.3133^ and 
K = 19, 18, 17, • • • to accuracy of two percent or less. 
Agreement becomes almost perfect if one notes that an 
0(E 2 ) correction makes the energy Ae loc of the localized 



state larger than A/(2ir£ 2 ) by factor / = 1 
includes this factor on the right-hand side of Eq. (3.2) 



\R 2 and 



The smooth interpolating current-field curve in Fig. 1 
rises from 10~ 6 to 1CP 2 as R varies from 0.23 to 0.33; this 
shows that derealization proceeds abruptly over a small 
range of R or E y . Actually, a simple parameterization 

tIoc 



J x oc = -e(E y /B)(l/L x )F n=0 (R) with 

F n=0 (R) = l/exp[exp[5.264- 11.52i?]] 



(3.4) 



gives a good fit to the curve for 0.17 < R < 0.33. 

Because the onset of the interference is determined ge- 
ometrically, the sequences R dlsc for analogous localized 
states of higher levels are essentially the same as the 
n = sequence, with slight deviations (of a few per- 
cent) in overall magnitude caused by the 0(E 2 ) correc- 
tion to Ae loc , which for the nth level is given by factor 
/„ = 1 + (n+ \)R 2 - The interference, however, becomes 
stronger with n (because of larger spatial overlap among 
nearly-degenerate states) and the interpolating current- 
field curves rise more steeply with R, as shown in Fig. 1 
for the n = 1,2,3 levels. These curves are distinct, but 
almost coincide by simple rescaling of the argument: 

F n=0 (R) & F n=l (R/ 1.07) 

« F n=2 (R/l.ll) « F n=3 (iJ/1.15), (3.5) 

which holds quite well for 0.001 < F n <• 0.01. One may 
define a critical ratio R cl as such at whi ch F n (R) reaches 
some given value, say, 1/10 3 . Then, Eq. (3.5) implies that 
R cr decreases only by a few percent as one moves to the 
adjacent higher level; indeed, from Fig. 1 one can read 
off R cr « (0.29, 0.27, 0.26, 0.25) for n = (0, 1, 2, 3). This 
shows that R — e£\E y \/(suo) is a good measure to express 
field-induced derealization of electron states, with only 
slight reference to specific levels n. 

We have made the sample twice and four times longer 
(L x — > 2L X — > &L X ) and observed that F n (R) scale with 
L x in such a way that (l/L x )F n (R) stay fixed for each 
n = 0, 1, 2, 3 over the range 0.17 < R < 0.32, giving rise 
to a relation analogous to Eq. ( |3.5| ) . It thus appears that 
F n= o(R) has a universal meaning of its own. 



IV. CASE OF RANDOM IMPURITIES 

Let us next consider samples with random impuri- 
ties. We earlier considerecta 180 short-range impurities 
of varying strength within the range \si\ < 0.1, randomly 
distributed over a domain of length ^/2n£ x 28 w 70^ 
and width \/2~k£ x 6 rs lh£. Here we shall reshuffle them 
and generate a new configuration of 360 impurities dis- 
tributed over a domain of larger width rs 2>0£\ see the ap- 
pendix for details. We embed this impurity configuration 
in the region — Ih £ < y < 15 £ oi a. Hall sample of length 
L x = J2ml x 28 « 70£ and width L y = V^£ x 18 « 45^, 
supporting N s = 28 x 18 = 504 electron states per sub- 
band with —\L y < yo < \L y . Note that the disordered 



domain accommodates roughly N° s rs 28 x 12 = 336 
electrons per subband. 

Of our main concern are the electron states residing 
over the disordered domain, that are to simulate electrons 
in the bulk of a realistic sample. We take no explicit ac- 
count of the edge states and simply leave the two zones 
15 £ < \y\ < 22.5 ^ of the sample impurity- free. In prac- 
tice, to determine the electron states in the disordered 
bulk one has to numerically diagonalize the Hamiltonian 
for a sample of somewhat larger width like this. 

For each subband n we arrange the electron states in 
the order of descending energy and use the number of 
vacant states, N v , to specify the filling of the subband. 
For convenience we start with the filled subband (N v = 0) 
and study the total current this subband supports, 



J, 



= I dyJ2i a) 



(x,y), 



(4.1) 



by removing electrons one by one. In Fig. 2(a) we plot 
J x as a function of N v for the n = subband and R = 
10~ 4 , 0.001, 0.01, 0.025, 0.05, 0.1, 0.2, where 



R= -ei E y / (0.1 w) 



(4.2) 



is now normalized using |s| max = Max{|s.;|} = 0.1 . 
There, simply for ease of exposing numerical results we 
have focused on 437 electron states whose center po- 
sitions (a\y\a) lie across the region \y\ < 19.6 £; the 
electron states residing outside are close to undisturbed 
modes. 

Let us first look at the R — 10~ 4 case, where the up- 
per and lower plateaus are clearly seen for < N v < 170 
and 270 < N v , respectively, supporting about 330 local- 
ized states in total. There the linear decrease of J x for 
170 < iV v < 210 and 230 < N v < 270 is due to grad- 
ual evacuation of extended states from the impurity-free 
zones \y\ > 15^. A steep fall of J x around N v ~ 220 in- 
dicates that a large portion of the Hall current is carried 
by a small number of states residing in the sample bulk. 

The plateaus shrink rapidly with increasing \E y \, and 
become hardly visible for R = 0.2 in Fig. 2(a). Somewhat 
unexpectedly, the way the upper and lower plateaus di- 
minish with increasing R is virtually the same for higher 
n = 1,2,3 subbands as well, as shown in Fig. 2(b) for 
the n = 2 subband. Numerically, the number of states 
forming the plateaus, iV pl , decreases exponentially with 
R, with little dependence on n = 0, 1,2,3, as seen from 
Fis_2(c); iV pl « 330 e ~ R / 0M1 gives a good numerical 

mis 

As seen from Fig. 2(a), a filled subband carries the 
same amount of Hall current as in the impurity-free case. 
One may thus simply look into the uppermost occupied 
subband in studying the dependence of the Hall current 
on the Hall field. 

In Fig. 2(b) the plateaus are somewhat blurred when 
E y is very weak, R < 10~ 3 . This blur or instability is due 
to the interference among nearly degenerate states and 



grows prominent in higher subbands, especially in the 
regime of small R, where J x itself gets small. Such blur 
makes A pl underestimate the actual plateau width, and 
this explains why apparent deviations from the universal 
A pl — R characteristic develop for small R in Fig. 2(c). 

We have so far supposed changing the filling fraction 
v = N e /N B of Landau subbands by varying electron pop- 
ulation N e with B kept fixed. Let us next consider what 
would happen if one reaches higher subbands by reducing 
B while keeping N e fixed, as commonly done in exper- 
iments. Unfortunately it is impractical to change the 
length unit I — \j\feB continuously in numerical simu- 
lations. We shall therefore take an indirect route and ex- 
amine what happens if the magnetic field is changed from 
Bq to Bq/k. We suppose that the sample sizcE3 is defined 
in units of £q — \/\feBo of the original field Bq. Then, as 
B is varied, the disordered domain accommodates about 
A| ff ss 336/k states per subband so that A; mp /A° ff S3 k. 
Thus, reducing the magnetic field Bq — > Bq/k not only 
attains a larger filling factor v — * k v but also effectively 
makes impurities k times denser. 

Some remarks are in order here. We have so far 
counted the number of localized states through the 
plateau width A pl , which best visualizes how the Hall 
plateaus shrink with R. This A pl , however, tends to over- 
look (a small portion of) localized states hidden off the 
plateaus and becomes less reliable for large R, where iV pl 
gets small. Fortunately there is another measure, less di- 
rect but better suited for detecting the disappearance of 
Hall plateaus: One can directly count the number 7V loc of 
states that carry an almost vanishing amount of current 
per state, e.g., less than 1% of the unit (~e/L x )E y /B. 
(We adopt this choice of 1% below.) 

Figure 3(a) shows A loc of the uppermost subband, 
plotted as a function of R oc e£\E y \/oj, for various choices 
of B = Bq/k with k — 0.33 ~ 4. Actually the data refer 
specifically to the n = subband but, in view of what 
we have learned from the fixed- B cases, we suppose that 
they apply to higher subbands as well. One may thus 
regard, e.g., the k = 2 curve as referring to the n = 1 
subband (or the spin-split n = subband if the elec- 
tron spin is taken into account) with 1 < v < 2, which 
is reached from the k = \ curve referring to the n = 
subband with < v < 1 by reducing B = B — > jBq. 

Figure 3(b) tells us that, in the regime of dense impu- 
rities Aj mp /A| ff S3 k > 1, the A loc — R curves appear 
quite universal, exhibiting, like A pl , an exponential fall- 
off over a wide range R > 0.1, with slight dependence 
on k for its range of change 1 < k < 4. The apparent 
deviations, growing with k for R < 0.1, are due to the 
interference among nearly degenerate states. Thus, an 
approximate scaling law of the form 



JVoe~ 



(4.3) 



with N ss N h 



360, determines how A depends on 



E y and B over a wide range R > 0.1. 

As for the normalization Nq it is clear that Ao oc Nj 
in general, and that Ao 



in which case all impurities are expected to capture 
electrons in the E y — > limit. It is remarkable that 
Ao ss Ai m p holds well even though the ratio Ai mp /A° 
far exceeds unity in Fig. 3(b). 

The situation somewhat changes in the regime of dilute 
impurities, A imp /A s cff < 1. See Fig. 3(a) and (c). There, 
as impurities become dilute (k = 1 — ► 0.33), A loc rises 
more rapidly with reducing R and tends to saturate for 
R ~ 0; this is because all the impurities would capture 
electrons before R reaches zero. Thus, for dilute impuri- 
ties A imp /A s cff < 1, the A loc - R characteristics signifi- 
cantly deviate from the scaling form ( |4.3| ). Still the values 
of A loc in the near-breakdown regime (R, = 0.25 ~ 0.4) 
stay roughly the same, and actually are consistent with 
those of the A imp /A s cff > 1 case. 

A physical picture suggested by this stability of A loc 
for large R is that the electron states remaining localized 
in the near-breakdown regime are always governed by 
the same set of impurities of large strength |s,| ~ |s| max ; 
that is, one always encounters the same set of localized 
states in the near-breakdown regime. If this is the case, 
weaker impurities (|sj| <C |s| max ) would be irrelevant to 
the onset of the breakdown. 

To test this picture we have examined the case of a 
truncated impurity configuration where only 180 impu- 
rities of strength 0.05 < |s,| < 0.1 are retained, with the 
weaker ones suppressed. The resulting A loc — R charac- 
teristics, plotted with unfilled symbols in Fig. 3(d) for 
K = 1,2,3,4, are indeed consistent with those of the 
original A S3 A imp = 360 case for R > 0.25. Fur- 
ther truncation of impurities into an even smaller range 
0.075 < | a,- 1 < 0.1 entails no essential change. This sup- 
ports our picture and implies that for a general impu- 
rity configuration Ao is not simply given by Aj mp but by 
the effective number of impurities extrapolated from the 
density of strong impurities (|sj| ~ |s| max )- We have also 
examined a number of cases with Ai mp = 180 and 720, 
and arrived at essentially the same conclusion. 

A practical way to define the disappearance of a Hall 
plateau is to refer to a critical value of E y at which A loc 
or the ratio N loc /N e reaches a certain prescribed number. 
Then, it follows from the universal N loc ~R characteristic 
in the near-breakdown regime that the critical field scales 
like el |^ r |/(0.1 w) = C = const., i.e., 



z\Ey\ - 



_Clu/£(xB 3 / 2 , 



(4.4) 



A imp when 



imp 



A^ 



imp 



N c 



where we have isolated the dependence of E" on the 
sample-specific impurity strength |s| max - It is understood 
that |s| max now refers to a typical strength of (strong) 
impurities relevant in the near-breakdown regime. EH The 
scaling law |-E"| ex B 3 / 2 applies to each spin-split sub- 
band reached by varying B. 

To get an idea of the magnitude \E y T \ let us look at 
Fig. 3(b) and choose C S3 0.25, for which we count only 
10 localized states in the uppermost subband. This leads 
to a critical field of magnitude \E y r \ ~ 250 V/cm at the 
v = 1 plateau for typical values w ss 10 meV and (. ~ 



100 A, with |s| max = 0.1. Critical fields of this order_of 
magnitude appear to be consistent with experiments. lj 
There are some testable consequences of the scaling 
law. First it is of interest to study how the Hall-plateau 
width depends on v and E y . Let B^ =1 ^ denote the value 
of B at which the lowest subband is filled up, so that one 
can write the filling factor as v = B ( - v=1 ^ jB. Let us 
suppose that each subband contains the same number of 
localized states ~ \N Xoc in its upper and lower halves 
of the spectrum. Then a plateau develops near integer 
filling v ~ n if \nN s - N e \ < ±A loc is satisfied; i.e., the 
plateau lies over the range f_ < v < v + with v± deter- 
mined from n/v± = 1 T \N Xoc {v±]E y )/N e , where A loc 
depends on v and E y . When E y is very weak and local- 
ization is almost complete, i.e., N loc — » N s , this leads to 
v± w n ± ^ so that the plateau width 

l\v=v+-v-K,\ (4.5) 

at each plateau, as observed prperimentallyE3 at low tem- 
peratures. [The observation^ 2 ] of sharp steps connecting 
the plateaus would suggest that the case of dense impu- 
rities Ni mp /N a > 1 applies to realistic samples.] 

In contrast, in the near-breakdown regime N loc /N e -C 
1 (or for large E y ) we find from Eq. (4.3) the plateau 
width AB/B("=V = l/z/+ - l/i/_ of the form 



AB 



-B^ =1 \N /N e )e~^ 3/2E y 



(4.6) 



around v ~ n, where £ = (C/0.058)/E y 



,cr (z^~l) 



refers 
to the critical field E y r ^ v ~ ' at the v = 1 plateau. 
Thus, near the breakdown the width AB scales like 
(l/n)e~£ n Ey at v ~ n plateaus. A close look into how 
the width AB decreases with E y at different plateaus 
will reveal another manifestation of the E" ex B 3 ! 2 law. 
The number of visible plateaus is readily determined 
also. Suppose that the total current J x is kept fixed 
while B is varied, as normally done in experiments. Let 
jcr (, ~ ,-j en0 ^ e thg cr itical current at which the v ~ 1 
plateau disappears. Then the critical current at the v ~ 
n plateau scales with n like J" (!y ~ n - ) = rT x l 2 j" '. 
As a result, for given J x one will observe n max plateaus, 
while reducing B, where n max denotes the largest integer 
satisfying the equation 



n < (J« ^/J x f. 



(4.7) 



states. This geometrical nature of the onset and sub- 
sequent abrupt growth make the derealization process, 
when expressed in terms of e£E y /B, practically indepen- 
dent of the Landau-subband index n. Such n indepen- 
dence is largely carried over to the case of random impu- 
rities. We have seen, in particular, that in the regime of 
dense impurities Ni mp /N s > 1 (which presumably applies 
to realistic samples) the number of localized states, de- 
creasing exponentially with \E y \, obeys an approximate 
scaling law written in terms of e£E y /uj ex E y /B 3 ' 2 . De- 
viations from the scaling form arise in the regime of di- 
lute impurities. Nevertheless the results of our numerical 
analysis are neatly summarized in a physical picture that 
the breakdown is controlled by a fixed set of impurities of 
large strength |sj| ~ |s| max , that support localized states 
in the near-breakdown regime. This picture offers a sim- 
ple explanation for the empirical E y T ex B 3 / 2 law of the 
critical breakdown field, and leads to some testable con- 
sequences concerning the number and widths of visible 
Hall plateaus- It is also consistent with the observation 
of Cage et alB that the breakdown of the dissipationless 
current is spatially inhomogeneous. 

Experimentally the E y r ex B 3 / 2 law is established not 
only for each specific sample bu)_also for samples differing 
in carrier density and mobility.til This experimental fact 
indicates that the critical field E" is practically indepen- 
dent of both electron concentration ex N e and impurity 
concentration ex A/imp- This is consistent with the above- 



V. SUMMARY AND DISCUSSION 



mentioned physical picture or its outcome in Eq. (4.4), 
where the critical field, though sensitive to the charac- 
teristic strength |s| max of disorder, appears practically 
independent of N e and Ai mp . The |s| max is related to the 
width of each subband. Its magnitude |s| max ~ O(0. 1) 
differs from a sample to another, but its range of varia- 
tion is presumably of 0(1), far smaller than the ranges of 
variations of sample-specific electron and impurity con- 
centrations. The observed E" ex B 3 / 2 law would thus 
be attributed to relatively small variations in |s| max for a 
variety of samples. It is also a natural consequence of the 
present intra subband process that the magnitude of E y T 
falls within the observed range of a few hundred V/cm, 
one order-of-magnitude smaller than what intersubband 
processesli-30 typically predict. 

We have used a uniform field E y to detect the current- 
carrying characteristics of Hall electrons. In real sam- 
ples the Hall field generated by an injected current is not 
necessarily uniform and has actually been observed to be 
stronger near the sample edges.E 3 ] Under such general cir- 
cumstances the field E y employed in our analysis should 
refer to the smallest of local averages of the Hall field in 
the sample bulk. 



In this paper we have studied the competition of the 
Hall field with disorder in triggering the breakdown of 
the QHE. In the case of a single impurity we have seen 
that field- induced derealization of a localized state starts 
with its interference with nearly-degenerate extended 
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APPENDIX A: IMPURITY DISTRIBUTION 

In this appendix we describe how to construct the 
configuration of 360 random impurities used in Sec. IV. 
We earlier used a configuration Pigo = {(x~i,yi,Si),i — 
1,---,180} of 180 impurities of strength Sj located at 
{xi, Vi) randomly distributed over the range < x < 24, 
< y < 6 and |s,| < 0.1. Generate out of this Pigo an- 
other set P'iso = {(xi + i 8 o, yi+iso, Si+i8o)j * = 1, ■ ■ ■ ) !80} 
via a shift (10,2,0.08) [which is our arbitrary choice] 
and rearrangement so that Xi+iso — x~i + 10 mod 24, 
l/i+180 = j/181-i + 2 mod 6, and s i+ i 8 o = Si + 0.08 mod 
0.2, within the range 0<x<24, 0<y<6 and \si\ < 0.1. 

Then shift Piso by (0, - 6, 0) and combine it with P'iso 
to form a set P360 of 360 impurities distributed over the 
wider range < x < 24 and — 6 < y < 6. Finally, via 
rescaling X{ = \/2tt£ x (28/24) x Xi and yt = ^/2~n£ x y~i 
we distribute these 360 impurities at positions (xi,yi) 
over the domain < x < 70^ and —151 < y < 15£. 
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FIG. 1. A potentially localized state of the n = 
level carries an amount of current J x oc [plotted in units of 
— {e/L x )Ey/E\, that becomes sizable at some discrete values 
of R — e£\Ey/siu\, with a smooth curve interpolating through 
the local minima. Also shown are analogous interpolating 
current-field curves for the localized states of the n — 1,2,3 
levels. 
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FIG. 3. The number of localized states per subband, iV loc , 
plotted as a function of R — el |i?j,|/|s|max us. 

(a) iV loc vs R with k = 0.33 ~ 4; 10 x N loc shown for R > 0.1. 

(b) logAf loc vs i?with k= 1,1.5,2,3,4. 

(c) k= 1,0.67,0.5,0.33. 

(d) log N loc vs R for a truncated configuration with only 180 
impurities of srtength 0.05 < s;| < 0.1 kept. 



FIG. 2. (a) The Hall current J x carried by the n — 
subband [in units of ~(e/L x )E y /B] vs vacancy N v for 
R=e£ |JS„|/|s| ma xW = 0.0001, 0.001, 0.01, 0.025, 0.05, 0.1, 0.2 . 

(b) J x vs N v for the n = 2 subband. The blur about the 
plateaus shown is associated with the R — 0.001 case, and 
the R — 10~ 4 curve is only partially shown. 

(c) The plateau width N pl decreases exponentially with R, as 
shown in a log iV pl vs R plot for the n = 0, 1, 2, 3 subbands. 
The lines joining the data points are meant to guide the eye. 



